Butterfly hysteresis loop and dissipative spin reversal in the S=l/2, V15 molecular 

complex 
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Time resolved magnetization measurements have been per- 
formed on a spin 1/2 molecular complex, so called V15. De- 
spite the absence of a barrier, magnetic hysteresis is observed 
over a timescale of several seconds. A detailed analysis in 
terms of a dissipative two level model is given, in which fluc- 
tuations and splittings are of same energy. Spin-phonon cou- 
pling leads to long relaxation times and to a particular "but- 
terfly" hysteresis loop. 

PACS numbers: 75.50.Xx, 75.45. +j, 71.70.-d 



In this letter we study the dynamics of the magnetiza- 
tion reversal of a molecular crystal made of nanometric 
molecules with non-interacting S — 1/2 spins. Despite 
the absence of energy barrier against spin reversal, this 
system shows hysteresis. This result are interpreted in 
details assuming spin rotation in a phonon bath, which is 
different from the situation of large spin molecules where 
only the spin bath is believed to be relevant 0-|^]. Res- 
onant phonon transitions are irrelevant, unless between 
states at different energies [pf or in the presence of a 
transverse field large enough to create a tunnel splitting 
of the order of the temperature energy scale o . 

The molecular complex K 6 [Vf^As 6 042(H20)]-8H 2 
(so-called V15) M is made of molecules with fifteen V IV 
ions of spin S = 1/2, placed in a quasi-spherical lay- 
ered structure formed of a triangle, sandwiched by two 
hexagons. The symmetry is trigonal (space group i?3c, 
a = 14.029 A, a = 79.26°, V = 2632 A 3 ). The unit-cell 
contains two V15 clusters and it is large enough that dipo- 
lar interactions between different molecules are negligible 
(a few mK) . All intra-molecular exchange interactions be- 
ing antiferromagnetic, the total spin of this molecule is 
S = 1/2. Such a small spin has zero energy barrier and 
relatively large splitting in zero applied field (~ 10~ 2 K). 
Although spin entanglement results in 2 15 eigenstates per 
molecule, the magnetization curves will be interpreted in 
terms of a dissipative two level model |q--|lO|] . 

Time-resolved magnetization measurements were per- 
formed with the micro-SQUID technique (50 — 400 mK, 
— 0.7 T/s) p!l| . In order to maximaze thermal con- 
tact with the bath, we choose a sample holder made by 
greece and silver powder and a small crystal of the V15 
(~ 50 /im). As an example we give a few hysteresis loops 
in Fig. [Da and Fig. 0a (only the positive parts are rep- 
resented, the other ones being rigorously symmetrical). 



When the field increases, coming from the negative satu- 
ration, the magnetization curve passes through the origin 
of the coordinates, reaches a plateau and then approaches 
saturation. This leads to a winged hysteresis loop char- 
acterized by the absence of irreversibility near zero field. 
Nevertheless, the initial susceptibilities being larger the 
faster sweeping field, the magnetization is out of equilib- 
rium also near zero field where it appears to be reversible. 
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FIG. 1. Measured (a—top) and calculated (b— bottom) hys- 
teresis loops for three temperatures and for a given field 
sweeping rate 0.14 T/s. The plateau is more pronounced at 
low T. The inset is a schematic representation of a two-level 
system Sz — ±1/2 with repulsion due to non-diagonal matrix 
elements. In a swept field the switching probability P is given 
by the Landau-Zener formulae (see text). The two levels are 
broadened by the hyperfine fields and the absorption or the 
emission of phonons can switch the polarization state of spins. 



The wings depend sensitivel 
field sweeping rate r. In Fig. [l 



on temperature T and 
where three hysteresis 



loops are presented at three different temperatures for a 
given sweeping rate, the plateau is higher and more pro- 
nounced at low temperature. The same tendency is ob- 
served at a given temperature and faster sweeping rates 
(Fig. 0a). When compared to its equilibrium limit (dot- 
ted curve in Fig. |2|), each magnetization curve shows a 
striking feature: the plateau intersects the equilibrium 
curve and the magnetization becomes smaller than at 
equilibrium. Equilibrium is then reached in higher fields 
near saturation. 
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FIG. 2. Measured (a— top) and calculated (b— bottom) hys- 
teresis loops for three field sweeping rates at T — 0.1 K. The 
observed plateau is more pronounced at high sweeping rate. 
The equilibrium curve can be approximated by the median 
of the two branches of the low sweeping rate hysteresis loop 
(dotted curve). In the top inset is plotted the spin and phonon 
temperature T s = T ph for T = 0.1 K and r = 0.14 T/s, when 
the field is swept from negative values. Ts decreases until 
zero-field and then increases linearly within the plateau re- 
gion. Then it overpasses the bath temperature to finally reach 
the equilibrium. In the bottom inset the calculated number of 
phonons with Ttw = Ah is plotted vs. the sweeping field mod- 
ulus (note the arrows) at equilibrium (T p h — Ts — T, dashed 
line) and out-of-equilibrium (tit h = nr=T s , r — 0.14 T/s, 
black line) . The difference between the two curves (thick seg- 
ment Aw) suggests the moving hole in the phonon distribu- 
tion, while their intersection gives the plateau intercept of the 
equilibrium magnetization curve. 

In order to interpret this magnetic behavior of the 



V15 molecules, we will analyse how the level occupa- 
tion numbers vary in this two level system (see Fig. [Lp 
inset) when sweeping an external field. In the ab- 
sence of dissipation, a 2-level model is well described 
by the bare Landau-Zener model, in the adiabatic or 
non-adiabatic case (low or high sweeping rates). The 
probability for the |l/2, -1/2) <-► |l/2, 1/2) transition is 
P = 1 — exp(— 7rAQ/4fi./iBr). In such a Landau-Zener 
transition, the plateaus of Fig. || should decrease if the 
sweeping rate increases, which is contrary to the experi- 
ments. Taking the typical value r = 0.1 T/s and the zero- 
field splitting A = 0.05 K |l^-jlq|, one gets a ground 
state switching probability very close to unity: in the 
absence of dissipation the spin 1/2 must adiabatically 
follow the field changes. Extremely large sweeping rates 
(f=a 10 9 T/s) would be needed to get into the quantum 
non-adiabatic regime P < 1. The mark of the V15 sys- 
tem is that the dissipative spin-phonon coupling is acting 
also near zero applied field because Tito rs Ao is of the or- 
der of the bath temperature, which is not the case for 
large spin molecules where Ao << fc^T. The spin tem- 
perature Ts is such that ni/712 = exp(A# /fee Ts), where 
Ah = vA~o + (2mb-Bo) 2 is the two levels field-dependent 
separation, and ni t 2(ni,2eq) the out of equilibrium (equi- 
librium) level occupation numbers. In the magnetization 
curves at 0.1 K (Fig. plpk), the spin temperature is sig- 
nificantly lower than the bath temperature T (n\ > rii eq , 
Ts < T) between —0.3 T (when the magnetization curve 
departs from the equilibrium one) and 0.15 T (the field 
at which the magnetization curve intersects the equilib- 
rium one). After this intersept Ts is larger than the bath 
temperature (ni < ni eq , Ts > T), and at sufficiently 
high fields (about 0.5 T) it reaches the equilibrium value 
(m = m^, Ts = T). 

In a direct process, the spins at the temperature 
Ts should relax to the phonons temperature within a 
timescale T\, the phonons being at the bath tempera- 
ture. However, even with a silver sample holder, it is 
not possible to maintain the phonon temperature equal 
to the temperature of the bath. This is because in V15 
below 0.5 K, the heat capacity of the phonons C p h is 
very much smaller than that of the spins Cg, so that 
the energy exchanged between spins and phonons will 
very rapidly adjust the phonons temperature T p h to the 
spin one Tg. Furtheremore, the energy is transfered from 
the spins only to those phonon modes with Hu> — Ah 
(within the resonance line width). The number of such 
lattice modes being much smaller than the number of 
spins, energy transfer between the phonons and the sam- 
ple holder must be very difficult, a phenomenon known 
as the phonon bottleneck jig]. Following |17| , the num- 
ber of phonons per molecule available for such resonant 
transitions is tit = /» cr(w)dw/(exp(?iw/fcT) — 1), where 
a(u))doj — 3Vlu 2 cIuj /(2ir 2 v 3 ) = number of phonon modes 
between u> and u> + du> per molecule of volume V, v is the 
phonon velocity and Aw is the transition linewidth due 
to fast hyperfine field fluctuations (they broden both en- 



ergy levels) fll8[] . Taking the typical values v » 3000 m/s, 
T w 10" x K and Aw « 5 • 10 2 MHz we find n T of the 
order of « 10~ 6 to 10~ 8 phonons/molecule. Such a small 
number of phonons is very rapidly absorbed, burning a 
hole of width Aui in the phonon density of states at the 
energy Tiuj = Ah [fL6|- If this phonon density of states 
does not equilibrate fast enough, the hole must persist 
and move with the sweeping field, leading to a phonon 
bottleneck. 

Now this description will be made quantitative. For 
a given splitting Ah, the time evolution of the two lev- 
els populations m,2 and of the phonon numbers n-r h at 
T p h obeys the set of two differential equations Of: (i) 
-hi = h 2 = P\ 2 n x - P 2 in 2 and (ii) h Tph = -(riT ph ~ 
n,T)/Tph — Pv2"n\ + P 2 \n 2 , where P\ 2 . 2 \ are the transi- 
tion probabilities between the two levels (they are them- 
selves linear functions of ut h ) and r p h ~ L/2v is the 
phonon-bath relaxation time (L is the sample size). Us- 
ing the notations x — (ri\ — n 2 )/(ni eq — n 2eq ),y — 
( n T ph ~ nT)/{riT + n/2) with n — J. <j(io)duj we get: 
(i) x = (1 — x — xy)JTi and (ii) y = —y/r^ + bx, where 
b = Cs/Cph and 1/ri = Pi 2 +P 2 i the direct spin-phonon 
relaxation time. By solving numerically this system for 
typical values, e.g. t\ = 10~ 2 s, r p h < 10 -6 s, b > 10 5 , 
we can see that T p h — » Ts 7^ T (phonon bottleneck) 
very rapidly, as expected. This leads to y — l/x — 1 and 
the second equation of the differential system becomes 
x = (x — £ 2 )/(l + bx 2 )/T p h- In the limit b >> 1 (in our 
case b rs 10 8 — 10 10 ) this equation has the solution: 



t/bTp h 



x + \n{{x-l)/{x -l)), (0.1) 



where xq = x(t = 0) and br p h is the spin-phonon recovery 
relaxation time (T p h = Ts — > T). When the system is 
not far from equilibrium (x ~ 1), we get an exponential 
decay of the magnetization, with the same time constant 
th = br p h- For a spin 1/2 system p7| : 



t h = a- 



tanh 2 (A H /2k B T) 



A 2 



(0.2) 



with a — 2n 2 h v 3 NT ph /3Au (N the molecule density). 

The dynamical magnetization curves calculated in this 
model are given Fig. plb and Fig. 0b. We started from 
equilibrium (xq = 1) in large negative fields. Then we 
let the system relax for a ver y short time St and we cal- 
culated x(5t) using Eq. [O]. This value was taken as 
the initial value for the next field (the field step is rSt). 
The parameters have been chosen to mimic the measured 
curves of Fig. |l]a and Fig. ga |19[. The obtained simi- 
larity supports the possibility of the phonon bottleneck 
effect at the timescale of a few 0.1 s. In the Fig. |p, in- 
set, we show the variation of the calculated spin-phonon 
temperature T s for T = 0.1 K and r = 0.14 T/s. We 
can note a linear variation in the plateau region (small 
positive fields, n\jn 2 w est.), after a cooling in nega- 
tive fields. The slope of this quasi-adiabatic linear region 
varies with the bath temperature and sweeping rate and 



gives the plateau dependence on these two parameters 
(see Figs. ||). 
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FIG. 3. The relaxation times th, measured (a— top) and 
calculated (b— bottom, same parameters as in Fig. hi 0b). 

In the Fig. gb inset we show the calculated field evo- 
lution of the number of phonons at energy Tiuj = Ah 
at equilibrium (T p h = Ts = T, dashed line) and out-of- 
equilibrium {jit h = nr=T s , r = 0-14 T/s, black line). 
The difference between the two curves (thick segment 
Alo) suggests the moving hole in the phonon distribu- 
tion, while their intersection gives the plateau intercept 
of the equilibrium magnetization curve (above which the 
hole dissapears and T p ^ — Ts > T). Let note that in zero 
field the system is out-of-equilibrium even if magnetiza- 
tion passes through the origin of coordinates (without a 
barrier, the switch between +1/2 and —1/2 follows the 
level structure shown Fig. [j] inset ). At larger fields, in 
the plateau region, Tu]n 2 w est. at timescales shorter 
than th = br p h (Eq. |0.2|) , even after the plateau crosses 
the equilibrium curve. Equilibrium is reached when th 
becomes small enough. 

Furthermore, we measured the relaxation of the mag- 
netization of our crystal at different fields and temper- 
atures, along the plateau region. The relaxation curves 
compared well to exponential decay and the obtained re- 
laxation times are presented in Fig. 0a. The comparison 
with those calculated (Fig. ||b) is acceptable. But we 



noted that a direct fit to Eq. 0.1 would necessitate larger 
values for a and A (« 0.4 - 0.6 sK 2 and w 0.2 - 0.3 K). 
Note that in V15 we have br p h > t\ and this leads to 
the phonon bottle-neck regime. However, in other sys- 
tems one might have br p h < T\ in which case the phonons 
would be at equilibrium but still with a butterfly hystere- 
sis loop (tjj is a linear combination of n and br p h p7[ ) . 
This type of hysteresis loop is general and characterizes 
dissipative spin reversal in the absence of barrier. 

In conclusion, the V15 molecular complex constitutes 
an example of dissipative two-levels system [B of meso- 
scopic size. The total spin 1/2 being formed of a large 
number of interacting spins, its splitting results from the 
structure itself of the molecule (intra-molecular hyper- 
fine and Dzyaloshinsky-Moriya couplings) and it is rather 
large (a fraction of Kelvin) p2]. In V15 and in other 
low-spin systems, splittings must be much larger than in 
large-spin molecules where the presence of energy barri- 
ers lowers them by orders of magnitude (e.g. 10~ n K 
in Mni2 JIB]). This is the reason why spin-phonon tran- 
sitions within the tuneling gap are important in low-spin 
molecules and not relevant in high-spin ones, unless a 
large transverse field is applied jq] (it increases the tun- 
nel splitting and probability) in which case we would also 
expect similar phenomena. 
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